Optical self-energy of superconducting Pb in the THz region 
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New THz data on the optical conductivity of Pb are presented as well as a detailed Eliashberg 
analysis with particular emphasis on phonon-assisted processes not included in a BCS approach. 
Consideration of the optical self-energy instead of the conductivity itself helps highlight the differ- 
ences with BCS predictions. Predicted coherence peaks are observed in the optical scattering rates. 
Impurities enhance the optical effective mass at zero frequency by an order of magnitude and induce 
a large peak at twice the gap in agreement with theory. This work illustrates the usefulness of the 
optical self-energy for the analysis of data. 



PACS numbers: 74.25.Gz,74.20.Fg,74.25.Nf,74.25.Kc 
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I. INTRODUCTION 

Optical properties of superconducting Pb have been 
measured in the past by several groups, including Palmer 
and Tinkham 1 -, who compared their data with BCS pre- 
dictions. Joyce and Richards^ saw phonon structure in 
the absorptivity and Farnworth and Timusk 3 analysed 
similar data to recover the electron-phonon spectral den- 
sity a 2 F(u>). This latter work was based on approxi- 
mate equations for the conductivity derived by Allen 4 
from Fermi golden rule considerations. Good agree- 
ment with the tunneling-derived electron-phonon spec- 
tral density was found 5 . More elaborate theories of the 
conductivity based on the Eliashberg equations and a 
Kubo formula for the current-current correlation func- 
tion were formulated by Nam^ and further elaborated 
upon in many places^ including Scalapino^ and Mar- 
siglio and Carbotte?' 1°. In this paper we report new 
high quality THz measurements on Pb. In addition to 
the new data, we provide a detailed theoretical analy- 
sis based on numerical solutions of the Eliashberg equa- 
tions with the known a F(u> ) ' . We emphasize par- 
ticularly the phonon-assisted processes not included in 
BCS and describe deviations that these processes intro- 
duce. BCS deals only with the coherent part of the sin- 
gle particle Green's function. In Pb, the electron-phonon 
mass renormalization factor A = 1.55. In this case, only 
1/(1 + A) = 1/2.55 of the spectral weight resides in the 
coherent part and A/(l + A) = 1.55/2.55 is found in the 
incoherent part which describes the phonon-assisted pro- 
cesses. 



II. METHODS 

The Pb films are deposited onto the (0001) surface 
of a single-crystalline sapphire substrate held at room 
temperature. Four-terminal resistance measurements on 
films of 30, 50, and 480 nm in thickness show the su- 
perconducting transitions around 7.2K with the width 



of ±50 mK, indicating little degradation in T c from the 
bulk. The present TDTS (Time Domain TeraHertz Spec- 
troscopy) measurements^ are carried out with the use 
of a commercial spectrometer (RT-20000, Tochigi Nikon 
Co. Ltd) based on standard transmission techniques, 
covering the frequency range of 0.15-2.0 THz. The real 
part <7i(w) and the imaginary part 02 (w) of the optical 
conductivity cr(uj) = ai(u))+i<j2(u>) are directly evaluated 
by numerically analyzing both the amplitude and phase 
of transmitted terawaves through the substrate (0.5 mm 
in thickness) with and without the film. The data will be 
presented for the 50 nm film most systematically studied. 

For a correlated electron system, specifically here a 
coupled, electron-phonon system, it has become standard 
practice to represent the optical conductivity a(T,u>) 
[with T, the temperature, and w, the photon energy] in 
terms of an optical self-energy in analogy to the quasi- 
particle self-energy which is introduced to include inter- 
actions in the one-particle Green's function. By defini- 
tion, 
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where lu p is the plasma frequency. The imaginary 
part of the optical self-energy gives the optical scat- 
tering rate 1/t°p(T,uj) = \T ', u>) and its real 
part gives the optical effective mass (m* p (T,uj)) accord- 
ing to w[m* p (T»/m - 1] = -2£° p (T,w). Further, 
m* p (T,Lo)/m — 1 is the optical mass renormalization 
A op (T, uS). Here, m is the bare electron mass. Tradi- 
tionally, it has been common to display in optical pa- 
pers the real and imaginary part of <j(T,u>) but in the 
recent literature, it has been recognized that much can 
be learned from the scattering rate and effective mass 
renormalization, particularly for comparison with quasi- 
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particle properties. In terms of <r(T, uj) 
1 ^2 / i \ 



T°P(uj) 



47T 



m*(T,uj) 



<r{T,u>) 



1 



47r erf 

, ,2 



(2) 



47r Im ( cr(T,w) 



(7 2 



4tt a? + ar^ 



The Drude conductivity is recovered when all interac- 
tions are neglected and only impurity scattering is in- 
cluded through a frequency and temperature indepen- 
dent scattering rate l/r™ p . In this case, the real part 
of £ op (T, uj) which must be related to its imaginary 
part by Kramers-Kronig is zero. For coupling to a sin- 
gle oscillator of energy in the normal state at zero 
temperaturi 
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where A is the coupling between electrons and phonons 
and 
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For uj — > 0, X op (uj) — 2A/lue = A, the quasiparticle 
electron-phonon mass renormalization. In this model, 
cr(uj) has two parts. One part is a delta function at 
uj = coming from the region uj < uje in ^ where 
there is no scattering. If impurity scattering is included, 
the (5-function broadens into a Lorentzian 
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FIG. 1: (color online) Real and imaginary part of the optical 
conductivity, <J\(T, u) and (72(T,u)), respectively, versus u) at 
various reduced temperatures t = T/T c . The points are the 
data and the lines are theory. 
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The second part which is finite only for uj > uje is a 
boson- assisted part: 
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Further, if we approximate in ([6]) m* p (uj)/m by 1 + A, its 
zero frequency limit which is exact for small uj, we obtain 
a coherent Drude part 
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with optical spectral weight of ujp/[8(l + A)] and effective 
scattering rate of 1 / [T lrnp ( 1 + A)] . The remaining spectral 
weight is to be found in the incoherent piece ([7]) and is 
uipX/[8(l + A)] which follows on application of the optical 
sum rule 
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(6) III. RESULTS 

In Fig. [H we present our experimental (points) and 
theoretical (lines) results for 6 temperatures as labelled 
by reduced temperature t = T/T c . In the top frame, we 
show the real part of the conductivity (T\(uj) in ri _1 cm _1 
versus uj up to a frequency of 2 x 10 12 Hz while in the 
lower frame we show the imaginary part a%(uj). For the 
real part, we note a kink in both data and theory at an in- 
termediate frequency which decreases in frequency, with 
increasing temperature. This arises because there are 
two separate absorption processes possible. At the low- 
est frequencies there can be direct absorption from the 
(8) thermally excited quasiparticles. This contribution pro- 
vides a Drude-like peak to o\(u>). At frequencies above 
twice the gap, direct absorption through the creation of 
two excitations out of the condenstate becomes possible. 
In addition, there is another incoherent phonon-assisted 
contribution to the absorption and this is the part of 
most interest in this work. However, it cannot easily be 
separated out of the real part of the conductivity. We 
will see that consideration of the optical self-energy S op 
(g) defined in Eq. |T]) helps in this regard. First, however, 
we note that the fit between theory and the data is very 
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good at low frequency in all cases. We have used the 
data in the normal state (t — 1.02) to determine the 
elastic scattering rate l/r lmp and the plasma frequency 
Lu p . Turning next to the imaginary part of a shown in the 
lower frame of Fig. [1] we note the same level of agreement 
and the points for t — 1.02 which are in the normal state 
also show good agreement with the theory. In general, 
the superconducting properties of Pb, while very differ- 
ent from those predicted by BCS, have been remarkably 
well-described by Eliashberg theory^ over a broad range 
of experiments and our new data shown here poses no 
exception. This is important given that in the past, it is 
has been argued that vertex corrections and other aspects 
of electromagnetic theory were needed to understand the 
data in Pb and yet this appears not to be the case here. 
With this level of agreement between theory and data, 
we should now be able to proceed to test for more sub- 
tle effects due to the phonon-assisted processes in the 
absorption. 

To see separately, more clearly, the role of the coher- 
ent absorption processes included in BCS theory and 
described by the coherent part of the electron Green's 
function G, and the incoherent phonon-assisted processes 
described by the incoherent part of G, it is convenient 
to deal with the real and imaginary parts of the opti- 
cal self-energy defined in Eq. {!]). Equations ([2]) and Q 
are for the optical scattering rate and effective mass, re- 
spectively. These are shown in Fig. [21 The top frame 
gives 1/t° p (uj) and the bottom [1 + X op (oj)] at three tem- 
peratures. To agree with the experimental data it was 
necessary to include an impurity scattering rate of 22 
meV. This is the only scattering that would be included 
in a BCS theory for which the normal state conductivity 
would be a pure Drude form. The optical scattering rate 
corresponding to this Drude is a constant flat line inde- 
pendent of frequency and temperature. We see clearly in 
the curve for the normal state that significant deviations 
from this constant occur. These are due to the boson- 
assisted processes which become even more important as 
the frequency is increased beyond the range of the data 
shown. This is shown in the top frame of Fig. [3l where 
we present our theoretical results over an extended fre- 
quency range up to 36 meV, well beyond the end of the 
phonon spectrum. For clarity, only the t = 0.52 (su- 
perconducting state) and t = 1.02 (normal state) curves 
are shown. Starting with the normal state (thick blue 
curve), we note a large deviation from the horizontal 
dashed black line at 22 meV for clastic impurity scatter- 
ing. The additional frequency-dependent scattering, due 
to boson-assisted processes, adds on to l/r mip according 
to Eq. ([?]) (in the case of a delta function spectral den- 
sity). For lj — > oo this contribution will saturate at 2ttA 
(see Eq. ©) even for a distributed a 2 F(uj) spectrum. 

Turning next to the superconducting state (thin red 
curve, top frame), we see larger boson structures with the 
vertical arrows indicating the transverse and longitudinal 
peaks in the phonon frequency distribution F(u>) of Pb 
shifted by 2A. These are preceded by a coherence peak 
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FIG. 2: (color online) The optical scattering rate l/r° p (T, ui) 
and mass renormalization [1 + X° P (T, ui)] vs ui. The points are 
experiment and the lines are theory. The data for t = 1.02 
are in the normal state just above T c and the other two data 
sets are below T c in the superconducting state. 



in the region above 2A which reflects the increased den- 
sity of final states for scattering in the superconducting 
as compared with normal state. This curve is strikingly 
different from the dot-dashed black curve which would 
apply to a BCS superconductor. An interesting compari- 
son is offered in the dotted curves which were also calcu- 
lated with the Pb a 2 F(u>) spectrum, but without residual 
scattering, i.e. l/r mip = 0. An important point to note 
is that no coherence peak is seen in the dotted thin red 
curve because this is an impurity effect, but at higher 
frequencies the dotted curve is very close to the solid one 
once it is displaced upward by 22 meV. This holds ex- 
actly for the normal state data, blue dotted (pure) and 
solid (with residual scattering). Returning to the lower 
frame of Fig. [2J we show the optical effective mass renor- 
malization [1 + X op (T, u)]. The solid blue dots apply to 
the normal state at t = 1.02 just above T c . Note that as 
uj — > A op ~ 1.7 which is just slightly larger than the 
zero temperature quasiparticle mass renormalization for 
Pb, which is 1.55. This is expected since it is known that 
the mass renormalization increases slightly with increas- 
ing temperature. Furthermore, the normal state effective 
mass has only a small frequency dependence predicted 
and this is observed. There is no dependence on impu- 
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FIG. 3: (color online) The calculated optical scattering rate 
1/t° p (T, ui) and optical mass renormalization [1 + \ op (T, io)] 
for a large energy range. Results are for our Pb calculation 
with impurity scattering rate of 22 meV. Red (thin line) is 
in the superconducting state and blue (thick line) is for the 
normal state, both at t = 0.52. We compare with BCS (dot- 
dashed and dashed black curves) and pure Pb (dotted curves) . 
The vertical arrows indicate the peaks in a 2 F(u) shifted by 
2A. 



rity scattering as it drops out of this curve. (This fact can 
be understood simply. The Kramers-Kronig (KK) trans- 
form of a constant is zero so that a constant scattering 
rate corresponds to no mass renormalization.) This is in 
striking contrast to the superconducting state for which 
impurities profoundly alter the optical curves. This is 
seen in the other two curves for t = 0.52 (red hexagons) 
and t — 0.82 (black squares). Note first the very large 
peaks which occur at twice the gap edge 2A. These peaks 
have not been emphasized before. They would not be 
there in the clean case (no impurities). Their origin is 
easily understood from our results for the scattering rate 
shown in the top frame. We see the sharp, almost verti- 
cal, rise in 1/t° p (uj) at 2 A with the energy scale of the rise 
given by the value of 1/ T imp = 22 meV. The KK trans- 
form of a step at is a logarithm with a singularity 
at wb. The larger the step, the larger the log singular- 
ity. Of course, here we do not quite have a vertical step 
and hence the corresponding peak in the effective mass 
is smeared. Another feature to note is the lj — > limit 



of X op . We note that Eq. Q can be rewritten as 



m* op (T,Lj) u 



LJ<T 2 



1)1 



4ir (coai) 2 + [wa^y 



(10) 



and that in the superconducting state in the limit of lj 
we get 
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where n is the electron density in the normal state and 
n s (T) is the superfluid density at temperature T. For 
T = 0, direct strong-coupling calculations of n s (T = 0) 
show that it is very nearly equal to n/(l + A) which leads 
to m* p (T = 0, lj = 0)/m ~ 1 + A, same as for the normal 
state, and reflects the fact that only the coherent part 
of the optical conductivity participates in the superfluid 
condensation. Returning to Eq. (fTU|) . one can verify that 
in the normal state, substitution of the Drude form © 
again gives m* op /m = 1 + A as we expect. If impuri- 
ties are included, as we have done in our work, n s is no 
longer n/(l + A) but in the dirty limit has instead the 
approximate form (in renormalized BCS)- 
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for a > 1 where a = l/[2Ar lmp (l + A)], which gives 
1 + X°p(T = 0, lj = 0) = 8.3. The exact numerical result 
(Fig. El bottom frame) is 9.5 at T/T c = 0.52. Most of 
the difference can be assigned to the small reduction in 
n s (T) between T = and t = T/T c = 0.52. 

In the lower frame of Fig. [3l we show additional theo- 
retical results on the optical effective mass which range 
over a larger frequency interval beyond that covered by 
the THz data. The dot-dashed black curve is the BCS re- 
sult (for l/r lmp = 22 meV) while the solid thin red curve 
includes as well the incoherent phonon-assisted Holstcin 
processes. We see clearly in this curve an image of the 
Pb phonon frequency distribution with transverse and 
longitudinal phonon peaks shifted by 2 A indicated by 
vertical arrows. The thin red dotted lines were obtained 
for pure Pb (i.e. l/r mp is set equal to zero and the 
only source of scattering comes from the inelastic part 
described by a 2 F(uj)). In this case there is no peak in 
[1 + X op ] at twice the gap but at higher frequencies the 
dotted and solid curves agree. For the normal state at 
t = 1.02, the solid thick blue curve applies. Note that for 
lj larger than a few times the phonon energy, normal and 
superconducting curves merge and all tend toward 1, i.e. 
X° p (lj) — > as lj — > oo. For BCS in the normal state, 
A op = for all lj as shown by the long-dashed curve. 



IV. A SUM RULE 

Recently, there have been several discussions of sum 
rules on 1/t° p in the literatur o 15 i 16 ' 17 . In particular, it 
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FIG. 4: (color online) The difference (l/rg P - vs w 

shown up to 35 meV and the partial sum, defined as I(ui) in 
the text. The black dashed curves are for the BCS case. In 
the inset, we compare data with theory. 



was noted that there exists a differential sum rule on the 
difference in optical scattering rate between the super- 
conducting and normal state. Defining the partial sum 
I{ui) as 



at t = 0.52 and shows a rapid increase through at twice 
the gap, after which a coherence peak is seen followed by 
boson structure beyond which the curve crosses the axis 
again. The solid thin red curve is the result for the partial 
sum I{oS). The dashed black curves are the BCS results 
for comparison. The crossing of the horizontal axis at a 
frequency slightly above the end of the phonon spectrum 
in the case for Pb (thin red curve) is due to the inelastic 
scattering processes and is not seen in the BCS curve. 
Note also the large positive peak. 

V. CONCLUSIONS 

We have presented new THz data on the optical con- 
ductivity of Pb and compared them with Eliashberg re- 
sults with an emphasis on phonon-assisted (Holstein) 
processes not part of BCS. We find that these show up 
most clearly in the optical self-energy rather than in the 
conductivity itself. In particular, the impurity compo- 
nent of the total scattering is seen directly in the op- 
tial scattering rate as a large constant which can easily 
be subtracted out to reveal directly the phonon-assisted 
part. In the normal state impurity scattering drops out of 
the optical effective mass leaving only phonon renormal- 
izations which show excellent agreement between theory 
and experiment. By contrast, in the superconducting 
state, the impurities provide a large renormalization of 
the zero frequency effective mass and cause a peak to ap- 
pear at twice the gap in the optical mass [1 + X op (T, u>)] 
which we predict here and observe in the experiment. 
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it can be showni& that I(fi — > oo) = for any value 
of \/T lmp and A. The energy scale for arriving at the 
asymptotic limit is set by A even in the dirty limit. As 
defined, in BCS-- Z(S7) is negative and has its maximum 
absolute value at 2A after which it begins to decrease to- 
wards zero, but always remaining negative. Our results 
for the case of Pb are shown in Fig. [4] and deviate signifi- 
cantly from the BCS behaviour just described. The solid 
thick blue curve gives the result for l/rg P (uj) — l/r^, p (w) 
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